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Abstract Phase dynamics has been measured in a driven mesoscopic Josephson oscillator 
where the resonance is tuned either by magnetic flux or by gate charge modula- 
tion of the Josephson inductance. Phenomena are analyzed in terms of a "phase 
particle picture", and by numerical circuit simulations. An analogy to switching 
of a DC-biased junction into voltage state is discussed. Operation principle of 
the recently demonstrated Inductive Single-Electron Transistor (L-SET) based 
on the driven oscillator is reviewed. The obtained charge sensitivity implies that 
a performance comparable to the best rf-SETs has already been reached with the 
novel device. 

1. INTRODUCTION 

Quantum measurement in the solid state has been shown to be feasible as 
several impressive realizations of qubits based on mesosocopic superconduct- 
ing tunnel junctions have emerged [1^1]. Sensitive measurement of physical 
quantities close to the limit set by the uncertainty principle is, on the other 
hand, an important issue in its own right. 

The Single-Electron Transistor (SET) is a basic mesoscopic detector, sen- 
sitive to electric charge on a gate capacitor. In order to gain advantage of the 
inherently large bandwidth (RsetCt,)^ 1 ~ 10 GHz of the SET charge de- 
tector, basically two new technologies have been developed where the SET 
is read using an LC oscillator built from "macroscopic" components but cou- 
pled directly to it. The "rf-SET" (Radio-Frequency SET) [5] is based on gate 
modulation of the Q-value of the oscillator. 

Because of limitations due to the dissipative nature of the rf-SET, a princi- 
pally non-dissipative "L-SET" (Inductive SET) technique has been developed 
very recently [6]. The L-SET is based on reactive readout of the Josephson 
inductance of a superconducting SET (SSET) in a resonator configuration [7]. 
Due to correlated Cooper pair tunneling, it does not exhibit shot noise or ex- 
cessive dissipation. 
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In this paper we first review the operating principle of the L-SET charge 
detector. Then we concentrate on discussing classical dynamics of the phase 
ip under the microwave drive, and we present new experimental data. In par- 
ticular, we discuss dynamical effects in the L-SET resonator which resemble 
the switching of a DC-biased Josephson junction into voltage state by drive or 
noise. 

2. THE L-SET CIRCUIT 

In the absence of DC bias voltage, a SSET has the Hamiltonian 
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(g ~ Igf 
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2Ej cos(<p/2) cos(0), 



(1) 



where q g is the gate charge, and ip is the phase difference across the SSET, as- 
sumed to be a classical variable here due to an environment having impedance 
much smaller than Rq ~ 6.5 kO. Ej is the single-junction Josephson energy, 
and the charging energy is related to sum capacitance by Eq = e 2 /(2(7s). 
Eigenvalues of this well-known Hamiltonian form bands E n (ip,q g ) (see Fig. 
!)• 




Figure 1. Two lowest bands n — and n = 1 from Eq. 1 for Ej/Ec = 1.7 (sample 2) 
plotted as function of the gate charge q 3 and the phase ip across the SSET. 

At the lowest band n = 0, the energy Eq increases approximately quadrat- 
ically as a function of the phase 99 when moved to either direction from the 
minimum at 99 = 0. This type of dependence is characteristic of an inductor. 
The effective Josephson inductance of the SSET is then 
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where the Josephson coupling has the effective value 

d 2 E (q g ,(f) 



(2) 



Ej 



dp 2 



(3) 



Here, $o = h/(2e) is the flux quantum. 

Since the energy band, and consequently, the inductance, depend on the gate 
charge q g , the L-SET electrometer is built so that the resonance frequency of an 
LC tank circuit connected to a SSET is tunable by q g . This allows in principle 
a purely reactive readout. 

The L-SET circuit we use is shown in Fig. 2 where the SSET is coupled in 
parallel to an LC oscillator resonant at the frequency fo = 1/(2tt)(LC)^ 1 ^ 2 , 
roughly at 600 MHz. The total system has the gate-dependent plasma reso- 
nance at f p = l/(2ir)(L tot C)^ 1 ^ 2 > fo, where L tot = L \\ Lj. The band- 
width Af ~ fp/Qe, where Q e is the coupled quality factor, is typically in the 
range of tens of MHz. 

The resistor r in series with the SSET is a model component for dissipation. 
As compared to the more standard way of drawing a resistor in parallel with 
the resonator, we got here a better agreement with the non-linear dynamics (see 
section 1.6.3). 



|r| V exp(io>t + arg(r)) 



Z, r 




Vo exp(icot) 



Figure 2. Schematics of the L-SET resonance circuit coupled to microwave feedline. The 
resonance is read by measuring amplitude or phase of the reflection coefficient T = (Z — 
Zo)/(Z + Zo), where Zo = 50Q, and Z(Lj) is the resonator impedance (including C c ). 

Charge detection is performed by measuring a change of amplitude or phase 
of the voltage reflection coefficient T, when the setup is irradiated by mi- 
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crowaves of frequency roughly f p . In the best sample so far [8], we have 
measured charge sensitivity 3 x 10~ 5 e/\/Hz over a bandwidth of about 100 
MHz. 

3. PLASMA OSCILLATIONS IN L-SET 

Dynamics of the L-SET oscillator can be analyzed in terms of a potential 
E p = E n (ip, q g ) + El due to the Josephson inductance and the shunting L, 
respectively, at the ground band approximately 

E p = -E (q g ) cos(<p) + $g/(8vr 2 L)(^ 2 . (4) 

At small oscillation amplitude, the phase particle experiences harmonic os- 
cillations around ip = 0, whose frequency f p is controlled by gate-tuning of the 
Josephson inductance (Fig. 3 (a)). This mode of operation, where the L-SET 
works as a charge-to-frequency converter, is the "harmonic mode". 




Figure 3. Illustration of the fictitious phase particle in a sinusoidally modulated parabolic 
potential. 

Once the oscillation amplitude is increased close to the critical current at 
ip = ±7r/2, the cos(ip) term changes the local curvature of the potential, and 
hence, the resonance frequency changes. At very high amplitude, Fig. 3 (b), 
the cosine wiggle becomes averaged out. Thus, we expect a change of resonant 
frequency from f p to /o roughly at an AC current of critical current magnitude. 
This change of resonance frequency when the sample is probed by critical 
power P c is reminiscent of a DC-biased Josephson junction switching into a 
voltage state [9]. 

At large excitations above P c , the highly nonlinear oscillator experiences 
complicated dynamics which does not in general allow analytical solutions. 
Numerical calculations over a large range of parameters, however, show con- 
sistently that the system response depends on L} also in this case [8]. This 
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mode of operation of the L-SET charge detector we call the "anharmonic" 
mode. 

4. SIMULATION SCHEME 

We simulated the transition from linear to non-linear oscillations in the L- 
SET circuit with the Aplac circuit simulation program which contains an im- 
plementation of a Josephson junction element. The SSET was modeled as a 
single tunable junction. We used the method of harmonic balance where am- 
plitudes of the first and three upper harmonics were optimized to create an 
approximate solution. Screenshot of the simulation schematics is shown in 
Fig. 4. 
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Figure 4. Screenshot from Aplac program used to simulate classical dynamics in the L-SET 
circuit (to aid convergence, we also used 1 kfl across the SSET). 



5. EXPERIMENT 

The measurements were done in a powerful Leiden Cryogenics MiniKelvin 
dilution refrigerator. At the base temperature of 20 mK, we used a microwave 
circulator which had about 20 dB backward isolation to cut back-action noise 
from the 4 Kelvin preamplifier [10]. 

The reflected probing excitation was amplified with a chain of amplifiers 
having a total of 5 K noise temperature, and detected with a network analyzer 
or a spectrum analyzer. 

We studied two samples (Table 1 and Fig. 5) where sample 1 had a tunable 
Ej. Since it had Eq — ksT, where the temperature was probably set by 
leakage through the circulator, its response had a hardly detectable gate modu- 
lation. Accordingly, we were able to study the plasma resonance in an almost 
classical junction. 

Sample 2 (discussed in Ref. [6], see section 1.6.1) was a sensitive detector, 
with a 15 MHz gate shift of the resonance frequency. 
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Table 1. Parameters of the two samples and of their tank LC oscillators as discussed in the 
text. Ej = hA/(8e 2 \Rset) is the single-junction Josephson energy, and Ec = e 2 /(2Ce) 
is the charging energy. For sample 1, SQUID structure of the individual SET junctions allowed 
tuning of the Ej /Ec ratio. 



sample 


Rset (kfi) 


Ej (K) 


E C (K) 


L*, (nH) 


L (nH) 


C*(pF) 


C c (pF) 


Q E 


1 


4.2 


3.5 ... 


0.17 


6 


3 


23 


0.72 


13 


2 


9.6 


1.6 


0.92 


16 


7.4 


8.4 


0.72 


18 




Figure 5. SEM micrographs of the samples discussed; sample 1 with SQUID junctions (left), 
and sample 2 (right) which had quasiparticle traps which were in contact to the electrodes 5/im 
from junctions. 



6. RESULTS AND DISCUSSION 
6.1 CHARGE DETECTION 

For sample 2, Ec 3> ksT, and hence, phase was properly localized in the 
potential. Accordingly, we got a rather good agreement for gate modulation 
of the resonance frequency (Fig. 6), where the resonance moves 15 MHz with 
respect to gate charge. This agrees well with the expectation of only 15 % 
modulation of Lj based on the rather high Ej/Ec — 1.7 of this sample. The 
minimum (w.r.t. gate) of Lj, fitted best with 18 nH is about 20 % higher than 
expected. This is probably due to a small amount of phase noise, or partially 
due to inaccuracies in determination of the parameters. 

Charge sensitivity was measured using amplitude readout. In the harmonic 
mode we got s q = 2.0 x lCr 3 e/\/Hz at the maximum power P c ~ — 116 dBm, 
corresponding to 20 fW dissipation in the whole resonator circuit. Note that 
due to Cooper pair tunneling, the power is not dissipated in the SSET island. In 
the anharmonic mode, significantly better sensitivity of s q = 1.4x 10~ 4 e/\/Hz 
was obtained. 
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Figure 6. Sample 2; measured frequency response at —125 dBm for successive gate values 
q g /e = 1 . . . 0, (a) amplitude, (b) phase. The insets portray the complete 2e periodicity [6]. 

In the best sample so far [8], having Ej/Eq — 0.3 we have measured 
3x 10 -5 e/\/Hz over a bandwidth of about 100 MHz. This kind of performance 
is comparable to the best rf-SET's [11], though power dissipation is more than 
two orders of magnitude lower. 



In Fig. 7 we plot data from sample 1 which behaves almost as a classical 
junction with negligible charging effect. By applied magnetic flux we tune 
simultaneously Ej of both SQUID-shaped junctions of the SSET, and hence, 
Lj. Periodicity of the resonance frequency with respect to applied magnetic 
flux is evident. 

A prominent feature in Fig. 7 is the absence of reflected signal at "critical" 
points where the resonance frequency switches to /o- This implies a strong 
dissipation which is visible as coupling to Zq. Pronounced dissipation at the 
critical points is representative of strong fluctuations at the artificial "phase 
transition" points. 

In Fig. 8 we plot f p measured at a low excitation, roughly 1/10 of critical 
current peak amplitude. The data are fitted to theoretical flux modulation 
of L*j. A symmetric SQUID has been assumed, where Ej tx cos(-7r<I>/<I>o) 
which affects the Ej/Ec ratio and thus Lj (see Eqs. 2 and 3). It is clear 
that something else is happening at low Ej, where the resonance meets /o 
more rapidly than expected. Rounding of the cusp in experimental data cannot 
be explained by asymmetry in the SQUIDS either, since then the resonance 
would deviate from true /o which is not the case. 

We argue this effect is due to "premature switching" caused by noise in the 
oscillator. The effect is again in analogy to physics in a DC-biased junction, 
namely, noise-induced switching and delocalization of phase. In contrast to 



6.2 



HARMONIC OSCILLATIONS 



8 




5.5 6 6.5 7 7.5 5.5 6 6.5 7 7.5 



f (Hz) x10 8 



Figure 7. Flux modulation of the power reflection (|F| 2 ) for sample 1 in order of increasing 
probing excitation (the graph labelled -103 dBm corresponds approximately to critical power 
P c at maximum Ej). 




Figure 8. Measured resonance frequency for sample 1 extracted from the first graph in Fig. 7 
(-123 dBm). The dotted line is a fit based on bare flux modulation of the Ej/Ec ratio. 



running of phase in a tilted washboard, the average of phase stays at zero in our 
case. Noise-induced switching happens when peak phase fluctuations reach 
(p c ~ tt/2. This condition corresponds to random motion of the phase particle 
at the bottom of the potential of Fig. 3 (a) with rms amplitude ~ (p c /2.5. 
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Noise excites the phase particle in the potential E p given by Eq. 4 up to an 
energy E p = ksT e jf, where T e ff is the electron temperature set by noise. 
Thus, approximately ksTejf = &l(ip c /2.5) 2 /(8Tr 2 L tot ). Since in Fig. 8 the 
switching happens at the flux ~ 0.8 x we get an estimate T e ff ~ 0.5 

K which is likely caused by leakage outside the band of the circulator. 

6.3 SWITCHING AND NONLINEAR OSCILLATIONS 

In this section, the discussion is based on experimental data from sample 
2. In Fig. 9, at P c ^ —116 dBm, the resonant frequency switches from the 
broad plasma resonance, centered at 720 MHz, into a narrower tank resonance 
at 613 MHz. The wavelike texture at —105 ... — 90 dBm is due to the cos(</?) 
Josephson-potential. 
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Figure 9. Measured frequency response (JF|, left) for sample 2 as a function of excitation 
power. "Critical power" P c (see text) is marked. Gate DC level was tuned approximately to 
the minimum of Lj, i.e., the maximum of f p . On the right is shown a corresponding Aplac 
simulation run using the circuit of Fig. 4. 

Changes in coupling, e.g., that Z goes through critical coupling at —105 
dBm, and sharpening of the resonance above P c , are due to an increase of the 
internal Q-value Qi from 20 up to several hundreds. Since Qi ~ uJoL*j/r, and 
since Lj acquires a large effective value at a high drive 3> 27r due to cancelling 
of positive and negative contributions, the supposed dissipation of r in series 
with the SSET has less effect at P > P c . 

The source of the dissipation modelled by the resistor r, probably situated 
in the SSET itself, is presently unknown, but it may be related to quasiparticle 
poisoning. It also limits the quality factor up to 20 approximately. 

Also plotted in Fig. 9 is an Aplac simulation for the circuit. It presents a 
qualitative agreement with experiment and predicts roughly correctly the end 
of the linear regime of plasma resonance. However, in experiment the switch- 
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ing is markedly sharp, which we attribute to the effect of higher bands of the 
SSET [12]. 
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